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This paper investigates the performance of a beta-configuration heat engine having a regenerative dis¬ 
placer. In the conventional beta-engine; the displacer and the power piston are incorporated in one 
cylinder. The displacer transfers the working fluid between expansion and compression spaces via the 
heater, the regenerator, and the cooler. In the present work, successive homogeneous layers of square 
wire meshes occupy the displacer space of a beta-engine that make the displacer to be a displacer and 
a regenerator simultaneously. The theoretical analysis of the engine is based mainly on Schmidt theory. 
The optimum dimensions of the heater, cooler, regenerator, piston stroke and displacer stroke as 
dimensionless ratios of the bore were found. The optimum phase angle between the piston and the 
displacer and the optimum ranges of the speed for each working gas were also found. In a comparison 
between the proposed engine which has a regenerative displacer and the GPU-3 engine which has 
a stationary regenerator and a solid displacer; it was found that; the proposed one delivers 20% more 
power with 10% more efficiency than the GPU-3 engine. 

© 2009 Elsevier Ltd. All rights reserved. 


1. Introduction 

The Stirling engine is an external combustion engine which uses 
compressible fluids as working media and is theoretically arranged 
to have Carnot efficiency. The thermal limits for its operation 
depend on the material used in its construction. In most instances, 
the engine operates at heater and cooler temperatures of 923 K and 
338 K respectively, [1]. Fig. 1 shows the p-V diagram of the ideal 
Stirling cycle with complete regeneration. The rejected heat during 
processes from state 4 to state 1 is absorbed during processes from 
state 2 to state 3. An external heat sink and heat source are required 
during processes 1-2 and 3-4 respectively, [2]. The first Stirling 
machine was built and tested by Stirling (1816) as an alternative to 
steam engines and other heat engines. Since Stirling engines are 
externally heated, they are powered using wide varieties of fuels 
and heat sources. The Stirling refrigerator does not have ozone 
depletion potential, [3]. Fig. 2 indicates three different configura¬ 
tions; namely alpha-, beta- and gamma-configurations. Each one 
has the same thermodynamic cycle, but it has different mechanical 
design characteristics. The alpha-configuration has twin power 
pistons in two cylinders and the working gas fluctuates via the 
heater, the regenerator and the cooler. In the beta-configuration; 
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the displacer and the power piston are incorporated in the same 
cylinder. The displacer transfers the working fluid between the hot 
space and the cold space through the heater, the regenerator, and 
the cooler. The gamma-configuration uses two separate cylinders; 
one for the displacer and the other is for the power piston. The 
power piston cylinder is connected to the displacer cylinder by 
using a connection port, [4]. The regenerator effectiveness depends 
mainly on its surface area, its void volume and flow rate through it, 
[5-7]. The laminate regenerator is fabricated from stainless steel 
screen sheets that were stacked on top of each other at certain 
angular orientation. The wire mesh regenerator is assembled from 
successive layers of stainless steel wire meshes of different pores/ 
inch. Heterogeneous arrangement of mesh layers enhances the 
regenerator effectiveness and reduces the pressure loss. A thermo 
hydraulic investigation on a packed bed solar heater having wire 
screen matrices of different geometrical parameters was mathe¬ 
matically modeled by [8]. The results indicated that the heater is 
thermo hydraulically efficient. The search for an engine cycle with 
high efficiency, multi-fuel and less pollution has led to reconsid¬ 
eration of Stirling cycle. Several engine prototypes were designed 
but their performances are still relatively weak if they were 
compared with those of combustion engines, [9]. This resulted from 
external and internal conduction, pressure drop in the regenerator 
and shuttle effect in the pistons. Current research indicates that 
solar powered, low temperature differential Stirling engines with 
gamma-configuration will be attractive in future, [10]. 
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Fig. 1 . Ideal Stirling cycle. 


This paper investigates the performance of beta-type Stirling 
engine using square wire meshes as regenerative media in the 
displacer space. The proposed engine with the regenerative dis¬ 
placer is shown in Fig. 3, where the working fluid fluctuates 
through a reciprocating regenerator instead of the stationary one. 
The selection of beta-type engine is referred to its higher specific 
power than alpha and gamma types, [11,12]. The selection of wire 
meshes as regenerative media is due to the rigidity under high 
pressure, high durability, stability under expansion and compres¬ 
sion and high regeneration effectiveness, [ 6 ]. The analysis of the 
engine depends mainly on Schmidt theory, [13], which deals with 
harmonic motion of the reciprocating elements, isothermal 
compression, isothermal expansion and perfect regeneration. 
The present study matched among the suitable dimensions of the 
heater, cooler and regenerator as dimensionless ratios of the 
cylinder bore. The piston stroke, the displacer stroke, the phase 
angle, the charging pressure and the speed of the engine were 
exhibited for six wire meshes. 


2. Beta heat engine 
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of fa bo is constant, i.e. 
r 3 / 2 ; hence 


0-1 = cos 1 [{r 3 /2 - (r g - r)cos fa/rfa 


(7) 


And as a consequence, the y-component of fa bo which is the 
displacer movement y b ' is: 


y b = V\ sin 0 i + r sin 0 


( 8 ) 


The proposed engine consists of single cylinder with hot head at 
its top, power piston, cooler and regenerative displacer as shown in 
Fig. 3. The proposed engine has porous regenerative displacer 
instead of solid displacer and fixed regenerator that are existed in 
the conventional beta-engine. Thus, the displacer performs two 
jobs; the displacement of the working fluid between expansion and 
compression spaces and the regenerator job. This was identified as 
an idea for gamma engine by [11 ]. The reciprocation of the pistons is 
controlled by using the rhombic drive mechanism shown in Fig. 4. 
The piston and the displacer are connected to cross bars c!b' and cb , 
respectively. The inter-meshing gears are used for the synchroni¬ 
zation of the reciprocation, [14]. The position and velocity vectors of 
each link of the drive mechanism are expressed as follows: 


Similarly, 

fa = cos - 1 [{r 3 /2 - (r g - r )cos <t>}/r 2 ] (9) 

The y-component of b , 0 which is the piston movement y b ,, is: 

= r 2 sin fa + r sin <f> ( 10 ) 

They-movements of both piston and displacer are represented in 
Fig. 5, where, one can get piston and displacer strokes and phase 
angle between them. Also, the swept volumes of both expansion and 
compression spaces can be stated. The variation of both expansion 
and compression spaces with the crank angle is expressed as 
follows: 
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Fig. 2. Different configurations of Stirling engines. After [4]. 
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Fig. 3. Scheme of the proposed beta-type Stirling engine. 


^ = (Vse/2) x (1+ cos 0) + V cU 


( 11 ) 


V c = (Vse/2) x (1 - cos 0) + (V sc /2) x [1 + cos(0 - a)] + V d>c + V P0 

( 12 ) 

The engine workspace is divided into three isothermal regions 
in the vision of Schmidt theory, [13]. The expansion and heater 
spaces are at high temperature T E , the compression and positive 
overlap spaces are at low temperature 7c and the regenerator space 
at an intermediate temperature T R which is expressed as follows: 


P = T 


m t RT E 


(Ye + V h ) + ll ]M(v R ) + l (Vc + V P0 ) 


(15) 


The Schmidt p-V ellipse is represented in Fig. 6. 

Referring to Fig. 3, the mass conservation equation for each 
space of the engine workspaces is as follows, [15]. 

For the heater space; 


dm h /dt 0 m h >E (16) 

For the expansion space; 


Tr = (T e - T c )/\n(T E /T c ) 


(13) 


According to Schmidt theory; the instantaneous pressure is kept 
constant throughout the engine workspace, and as a consequence, 
the mass of the gas is: 


m t 


p(V E + V h ) p(V R ) p(V c + Vpp) 
RT e RT r rt c 


(14) 


dm E /dt = m' h ^ E - m E _ R 
For the regenerator space; 

dm R /dt = m E _ R - m R ^ P0 

For the positive overlap and compression spaces; 

dm c /dt m R ^ P0 - 0 


(17) 

(18) 


The Schmidt pressure is: 


(19) 
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O Origin 

o' Crank shaft axis 

a Crank shaft pin 

ab Displacer connecting 

rod 

be Displacer cross bar 

ab' Piston connecting rod 

b'c' Piston cross bar 

r Crank shaft radius 

r Displacer connecting 

rod length 

r 2 Piston connecting rod y 

Length 4 

r 3 Cross bar length_ 

0 Crank angle 

(j\ Displacer connecting 

rod angle 

02 Piston connecting rod 

angle 

yb Displacer movement 

3V Piston movement 



Fig. 4. Dimensions and angles of rhombic drive. 


From equations (16)—(19); one can find; m* h _ >£ , rn E ^ R and 
m R ^ po . Hence the flow rate through the heater, the regenerator 
and the cooler can be found as follows: 

m h = m E ^ h ( 20 ) 


Re R = m R x <4yd, k/(Mk x q r) 

(24) 

Re k = m k x d k /(/i c x a k ) 

(25) 

2A. Heater 



( m E—>R + m R^ PC))/2 

(21) 

m R-+ PO 

(22) 


Fig. 7 shows the periodic flow rate through the heater, the cooler 
and the regenerator. The instantaneous Reynolds numbers through 
the heater, the regenerator and the cooler are found as follows: 

Re h = m h x H h /(n E x a h ) (23) 



The heater of the proposed engine is similar to the one which 
was experimentally tested by [16]. It is of inverted U-shaped cups 
that can be externally heated by flue gases as shown in Fig. 8. The 
experimental tests of such heaters concluded that the overall heat 
transfer coefficient during a complete cycle is; U = 311 W m 2 °C for 
exhaust of LPG fuel at temperature 7}= 600 °C, Thus; 


$ = u\ 0 (r f -T E ) 


(26) 


The pressure drop of the gas due to its flow through the heater 
is, [17]: 

Ap h = 1.6 x \Kxpxv 2 /2\, (27) 
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Fig. 5. Piston and displacer movement versus crank angle. 


Fig. 6. Schmidt p-V diagram at a charging pressure of 20 bar. 
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Fig. 7. Periodic mass flow rate through cooler, heater and regenerator. 

2.2. Regenerator 

The regenerative displacer has cylindrical shape, as shown in 
Fig. 8. It is formed of successive layers of stainless steel wire meshes 
having the same mesh size. The layers were homogenously stacked 
on the top of each other without revolving angle. Referring to Fig. 9, 
six meshes of different dimensions are suggested to be used in the 
proposed regenerative displacer; their technical data are given in 
Table 1, [6]. 

The pressure drop of the gas due to its flow through the dis¬ 
placer is, [18]: 


J R 


hpR = f X ( H/d hd ) X p X v 2 /2 
where 

<4yd = [WO - '/')] x d w 

The friction factor,/, is expressed as follows, [18]: 

log (f R ) = 1.73-0.93 log (Reg) 

If 0 < Re R < 60 

log (f R ) = 0.714 - 0.365 log(Ke R ) 

If 60 < Re R < 1000 

log (f R ) = 0.015-0.125 log (Reg) 

If Re R > 1000 

The regenerator effectiveness is, [18]: 
e = NTU r /(1 + NTUr) 

NTUr = 2 x St R x Hr/ d hyd r 
and 
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Fig. 8. Half cross-section in the heater and the cylinder head. 
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Fig. 9. Schematic diagram of wire mesh layer, After [6]. 


St r = 0.595/Reg 4 x Pr R 
2.3. Cooler 


(35) 


The water jacket about the compression space is the cooler in 
the proposed engine. The working fluid is cooled during its fluc¬ 
tuation through the compression space. The cooling rate is, [19]: 


Q/< — Mk,o(^o ^w) — htkkLk (Jk i /f< 7 o) 
= hiAj<j(T c - T ki ) 


(36) 


The pressure drop of the gas due to its flow through the cooler is, 
[17]: 


Ap fc = 1.6 (/L/d) xpxv 2 / 2 


Jk 


(37) 


(29) 3. Engine power and efficiency 


The elliptic p-V diagrams for expansion and compression spaces 
when taking the hydraulic losses into consideration are illustrated 
in Fig. 10. The indicated power is, [12]: 


Pind = [$Pe dVs + j>Pc 

The rate of heat added is: 

0 ' h = Nx$P E dV E 

The indicated efficiency is: 

Pind = ^md/Q/i 


(38) 


(39) 


(40) 


A computer program in the form of a spreadsheet was prepared 
on the Microsoft EXCEL to compute the instantaneous volumes of 
expansion and compression spaces, instantaneous mass flow rates 
through the three heat exchangers, instantaneous Reynolds 
numbers in the heat exchangers, heat transfer rates, pressure drop, 
indicated power and efficiency. This was done for a complete 
revolution of the crank at a step of 1.0°. It was assumed that the 


Table 1 

Technical data of wire meshes. 


Mesh-i 

D w (mm) 

Porosity; fp (%) 

Mesh-50 

0.200 

69.08 

Mesh-100 

0.112 

65.37 

Mesh-200 

0.050 

69.08 

Mesh-300 

0.030 

72.17 

Mesh-400 

0.030 

62.89 

Mesh-500 

0.025 

61.35 
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Fig. 10. Pressure-volume diagrams for both expansion and compression spaces. 
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Fig. 11. Variation of power versus heater height for different wire meshes. 
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Fig. 12. Variation of efficiency versus heater height for different wire meshes. 


engine works between temperature limits of r m ax = 600 °C and 
T min = 80 °C. The cylinder bore is 100 mm. The program has the 
possibility to vary the dimensions of the heat exchangers, piston 
stroke, displacer stroke, phase angle, charging pressure and speed 
for each mesh size. 


4. Results and discussion 

The heater height to bore ratio, Hh/D was changed from 0.75 to 
1.55 at a step of 0.1, the corresponding values of both power and 


Fig. 14. Variation of efficiency versus displacer height for different wire meshes. 



Sp/D 

Fig. 15. Variation of power versus piston stroke for different wire meshes. 


efficiency were determined. Fig. 11 shows higher rate of increase in 
the power with the increase of the heater height ratio up to about 
1.0. Beyond this value, the rate of increase of the power reduces 
some how. This is referred to the increase in the heat transfer area; 
however, further increase in the heater height causes an increase in 
the dead volume which reduces the power. Also, thin wire meshes 
cause high pressure losses than the thick ones that reduce the 
power as shown in Fig. 11. The same tendency is even remarked in 
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Fig. 16. Variation of efficiency versus piston stroke for different wire meshes. 


Fig. 19. Variation of power versus phase angle for different wire meshes. 



Fig. 17. Variation of power versus displacer stroke for different wire meshes. 



V, (m3) 


Fig. 18. Schmidt p-V diagrams at different displacer strokes for mesh 100. 

Fig. 12 for the efficiency. Therefore; the heater height was selected 
to be Hh/D = 1.0 for every mesh size. 

For the same displacer height, the thin wire meshes show more 
decrease in the power than that of thick ones as shown in Fig. 13. 


Table 2 

Optimum displacer stroke for each mesh size. 


Mesh size 
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Fig. 20. Variation of efficiency versus phase angle for different wire meshes. 

This is due to the increase in hydraulic losses. The increase in the 
displacer height increases the mass of the gas which increases 
the power. Also, the increase in the displacer height increases the 
regenerator effectiveness which increases the efficiency as shown 
in Fig. 14. 

Increasing piston stroke increases the mass of the charged gas, 
which increases the indicated power up to about S p /D = 0.5, as 
shown in Fig. 15. Beyond this value, the more increase in the piston 
stroke increases the charged gas which results in high pressure 
losses through the heat exchangers and as a consequence the 
power decreases. A reduction in the efficiency versus the piston 
stroke is noticed in Fig. 16, this is due to the increase in the heat 
rejection rate which resulted from the increase in the cooler surface 
area. Referring to Fig. 17, it is obvious that; each mesh size has an 
optimal value for the displacer stroke as given in Table 2. 

Increasing the displacer stroke increases the charged gas which 
increases the power. Further increase in the displacer stroke 
increases the hydraulic losses, moreover; more increase in 
hydraulic losses was noted for thin meshes. So, the optimum dis¬ 
placer stroke decreases as the mesh size increases. The increase in 
the elliptic p-V area due to the increase in the displacer stroke is 
clear in Fig. 18. 

The phase angle can be varied by the variation of the dimensions 
of the rhombic drive mechanism, [20]. The change in the phase 
angle causes a change of both compression and pressure ratios that 
affect the work per cycle and consequently the power. The phase 
angle was varied from 45° to 135°. Referring to Fig. 19, it is evident 
that; the more optimum phase angle for most of meshes is about 
90°. The optimum phase angle results in high temperature ratio, 
? which reduces the efficiency as shown in Fig. 20. 
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Fig. 21. Variation of power versus the speed for air as a working gas. 



N, (rpm) 

Fig. 22. Variation of efficiency versus the speed for air as a working gas. 


Fig. 24. Variation of efficiency versus the speed for nitrogen as a working gas. 
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Fig. 25. Variation of power versus the speed for hydrogen as a working gas. 
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Fig. 23. Variation of power versus the speed for nitrogen as a working gas. 


Figs. 21,23,25 and 27 show the effect of the speed on the power 
if the engine was charged with air, nitrogen, hydrogen and helium 
respectively. The effect of the speed was obtained using the 
optimum dimensions of heat exchangers, optimum strokes and 
optimum phase angle for each mesh size. The optimum range of the 
speed for each gas has an acceptable range of the efficiency as 
shown in Figs. 22,24,26 and 28. 

Regardless of the gas leakage and the lifetime of piston rings, 
one can say; the more charging pressure, the more power as shown 



N, (rpm) 


Fig. 26. Variation of efficiency versus the speed for hydrogen as a working gas. 

in Fig. 29. This is due to, the increase in the charged gas. The 
increase in the charging pressure causes a decrease in the regen¬ 
erator effectiveness as shown in Fig. 29. This will increase the 
temperature ratio, ? and as a consequence the efficiency decreases 
as shown in Fig. 30. 

5. Comparison between present engine and GPU-3 engine 

The urgent need to renewable energies has led to the develop¬ 
ment of Stirling engines that have an excellent theoretical effi¬ 
ciency, equivalent to Carnot. A previous modeling was applied on 
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Fig. 27. Variation of power versus the speed for helium as a working gas. 


Fig. 30. Variation of efficiency versus charging pressure for air as a working gas. 




Fig. 28. Variation of efficiency versus the speed for helium as a working gas. 


Fig. 31. Comparison between power per cc from present engine and GPU-3 engine. 



Pch ( bar ) 

Fig. 29. Variation of power versus charging pressure for air as a working gas. 


the General Motor GPU-3 Stirling engine by [9]. The GPU-3 engine 
has a swept volume of 233 cm 3 and uses helium at a mean pressure 
of 41.3 bar. The model considered T £ = 977 K and Tc = 288 K. The 
data of GPU-3 engine were fed to present model to calculate 
the corresponding results that were compared with those of [9]. 
The comparison shows that the present engine delivers about 20% 
more power per cc with about 10% more efficiency as shown in Figs. 
31 and 32. As the porosity increases the present engine develops 
more power at high efficiency than that of the GPU-3 engine. 



Fig. 32. Comparison between efficiency of present engine and GPU-3 engine. 

6. Conclusions 

The present work exhibited the performance of beta-type Stir¬ 
ling engine having regenerative displacer. Successive homogeneous 
stainless steel wire meshes filled the space of the displacer of 
engine. The porous displacer performs as a displacer and as 
a regenerator simultaneously. Six wire meshes of different mesh 
size were suggested to be the regenerative media. The analysis of 
the engine is mainly based on Schmidt theory. The optimum 
dimensions of heater, cooler, regenerator, piston stroke, displacer 
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stroke as dimensionless ratios of the cylinder bore were found. The 
optimum phase angle and the optimum ranges of the speed for 
each gas were also stated. The main concluded remarks of the 
present work can be systemized as follows: 

1. The proposed beta-engine which has a regenerative displacer 
instead of the solid displacer and the stationary regenerator 
delivers 20% more power per cc than that from the GPU-3 
engine which has a stationary regenerator and solid displacer 
with an increase in the efficiency of about 10%. 

2. The present model suggested successive layers of homogenous 
wire meshes of different mesh sizes. However, the model 
allows using of successive layers of heterogeneous meshes as 
well as other types of regenerative media. 

3. The replacement of the stationary regenerator by a recipro¬ 
cating regenerative displacer incorporated inside the engine 
cylinder of beta-type engines will reduce the space occupied by 
the engine and it will slightly reduce the engine weight. 

4. The present model exhibited the dimensions of the engine as 
dimensionless ratios of the bore. And as a consequence, the 
engine sizes can be adapted to be similar prototypes. 
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Appendix. Nomenclature 

A surface area, m 2 

a cross-section area, m 2 

cc cubic centimeter 

D cylinder bore, m 

d diameter, m 

f friction factor 

H height, m 

h heat transfer coefficient, W m -2 K -1 

i wire mesh, pores/inch 

i , j Cartesian unit vectors 

I< minor losses factor 

k thermal conductivity, Wm _1 K” 1 

L length, m 

m mass, kg 

m mass flow rate, kg s^ 1 

N speed, rps 

NTU number of transfer units 

P power, W 

p pressure, N nrT 2 

Q* heat transfer rate, W 

R specific gas constant, J kg K _1 

Re Reynolds number 

~r position vector, m 

S stroke, m 

St Stanton number 

T temperature, I< 

t time, s 

U overall heat transfer coefficient, W nrT 2 K _1 

V volume, m 3 

v velocity, m/s 

~v velocity vector, m/s 

w angular velocity, rad/s 

x,y cartesian coordinates, m 

a phase angle, rad 

s regenerator effectiveness 


</> crank angle, rad 

7] efficiency 

p viscosity, kg m s -1 

? temp, ratio = TdT E 

p density, kg m -3 

V/ porosity 

Subscripts 

C compression space 

cl clearance 

ch charging 

d displacer 

E expansion space 

f flue gases 

g inter-meshing gear 

h heater 

hyd hydraulic 

i inner conditions 

ind indicated 

k cooler 

o outer conditions 

p piston 

PO positive overlap 

R regenerator 

SE swept expansion 

SC swept compression 

t total 

w wire 
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